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We find the leakage current through a double quantum dot in the Pauli spin blockade regime 
accounting for inelastic (spin-flip) cotunneling processes. Taking the energy-dependence of this 
spin-flip mechanism into account allows for an accurate description of the current as a function 
of applied magnetic fields, gate voltages, and an inter-dot tunnel coupling. In the presence of an 
additional local dephasing process or nonuniform magnetic field, we obtain a simple closed-form 
analytical expression for the leakage current giving the full dependence on an applied magnetic field 
and energy detuning. This work is important for understanding the nature of leakage, especially 
in systems where other spin-flip mechanisms (due, e.g., to hyperfine coupling to nuclear spins or 
spin-orbit coupling) are weak, including silicon and carbon-nanotube or graphene quantum dots. 

PACS numbers: 72.25.Rb,73.63.Kv,73.23.Hk,85.35.Be 



I. INTRODUCTION 

Spin-dependent current blockade effects have been ob- 
served in quantum dots^ as well as in molecular— and 
single-atom tansport^ These effects are the basis of spin- 
to-chargc conversion schemes, essential for measurements 
of spin coherence and relaxation^— as well as the accurate 
initialization and readout of spin states for spin-based 
quantum information processing^ Blockade effects have 
further allowed the observation of intriguing slow peri- 
odic oscillations in current, dependent on nuclear spinsJ^ 
A detailed microscopic understanding of how this block- 
ade can be lifted is important to develop an accurate 
description of these effects and to point the way to gen- 
erate a more robust blockade for the study of further 
spin-dependent phenomena. 

The Pauli spin blockade of current through a double 
quantum dot occurs when each of two quantum dots 
in series energetically favors a one-electron configuration 
(we will refer to this as the (1,1) regime, where (n, m) 
refers to n electrons on the left dot and m electrons in 
the right). Restricting to only the lowest non-degenerate 
single-particle orbital state in each dot, there are four 
possible spin configurations in the (1,1) subspace: one 
spin singlet and three spin triplets. To generate sequen- 
tial transport of electron charge from left to right, the 
double-dot must pass through the (0, 2) charge config- 
uration, but due to the Pauli exclusion principle, the 
lowest-energy (0, 2)-statc is a spin singlet when only the 
lowest single-particle orbital state is accessible. An inter- 
dot tunnel coupling preserves the spin of the two-electron 
state and therefore couples only the (1, 1) singlet to the 
(0, 2) state. After a small number of electrons has passed 
through the double-dot, eventually one of the spin-triplet 
states will be occupied by chance, leaving the double 
dot stuck in a "blocked" configuration. This blockade 
can be lifted either through the direct hybridization of 



singlet and triplet states with the addition of spin-non- 
conserving terms to the Hamiltonian (due, e.g., to the 
spin-orbit or hyperfine interactions), or through direct 
energy-conserving transitions between triplet and singlet 
levels. In spite of this relatively simple explanation for 
the Pauli spin blockade, the situation is complicated by 
several possible microscopic mechanisms that may domi- 
nate in determining the leakage current depending on the 
material and device characteristics. It is therefore impor- 
tant to understand precisely what influence each of the 
possible microscopic mechanisms may have on the over- 
all leakage current in order to identify the most relevant 
mechanism and possibly to suppress it. 

Most mechanisms that lift spin blockade are particular 
to the materials used to manufacture a double-dot device; 
the contact hyperfine interaction between electron and 
nuclear spins lifts the blockade in GaAs double dotsr^ - — 
and a strong spin-orbit interaction plays the predomi- 
nant role in lifting the blockade in InAs nanowire double 
dotsJ^r— Both of these mechanisms can be suppressed by 
manufacturing double dots using siliconi^r— or carbon- 
base d 21 ' 22 materials, in which the majority isotope has no 
nuclear spin and the spin-orbit coupling strength cx Z 4 
is significantly weaker due to a smaller atomic number 
Z. One blockade-lifting mechanism that is present in all 
double-dot devices, independent of the material composi- 
tion, is exchange of spins with the leads through higher- 
order tunneling (cotunneling) processesi 17 : 23 " — By un- 
derstanding and controlling these processes, one can ac- 
curately calibrate single-spin readout and improve on 
rapid spin preparation schemes^ Moreover, cotunneling 
processes have been shown to be significant in determin- 
ing dynamic nuclear-spin polarizaiton processes, both in 
experiment^ and in theor y 26 : 27 so a further understand- 
ing of cotunneling may allow for the preparation of a 
more highly-polarized nuclear-spin system. 

In this article we derive analytic expressions for 
leakage-current lineshapes accounting for inelastic cotun- 
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neling processes in well-defined and generically accessi- 
ble limits. Some results of this analysis have recently 
been shown, experimentally, to be consistent with trans- 
port measurements on silicon double quantum dots^S 
and have proven useful in determining microscopic pa- 
rameters associated with those devices. A similar ap- 
plication of the results presented here to other material 
systems may shed light on, e.g., unusually broad line- 
shapes in the magnetic-field-depcndcnt current through 
isotopically enriched 13 C nanotube double dots^ 



Inelastic (spin-flip) cotunneling has been known as a 
significant spin-flip mechanism since early measurements 
of triplet-to-singlet decay in vertical double quantum 
dots,— where the triplct-to-singlet decay rate was shown 
to be limited by inelastic cotunneling^ 3 - In the context of 
the Pauli spin blockade, spin-flip cotunneling rates have 
been calculated and compared to experimental data in 
the high-temperature regime where the associated tran- 
sition rates between energy levels are independent of the 
energy-level spacing i 17 i 24 More recently, the consequence 
of the full energy dependence of these rates has been 
calculated 2 ^ and verified in experiment^ In this pa per 
we apply and extend the analysis presented in refs. [25| 
and [20| to a broader range of parameters and provide 
a general and intuitive formalism for the calculation of 
leakage current through blockaded structures. For sim- 
plicity, in specific calculations we neglect orbital/valley 
degeneracy in our treatment, which may be relevant for 
quantum dots made from graphene, carbon nanotubes, 
or silicon nanostructures and can lead, in general, to a 
more complicated spin-valley blockadei 29 i 30 However, the 
general formalism we present can also be applied directly 
to systems with valley degeneracy and many of the re- 
sults we present will be qualitatively unchanged in the 
presence of additional orbital degeneracies. 



The remainder of this article is organized as follows: In 
Sec. [II] we present an intuitive and general procedure for 
the calculation of current through blockaded structures 
given a set of decay rates obtained from a microscopic cal- 
culation. In Sec. IIIII we specialize to the case of a double 
quantum dot in the Pauli spin blockade regime. We re- 
call the calculation of sequential-tunneling and spin-flip 
cotunneling rates from ref. [25| and apply the procedure 
of Sec. |TT] to find leakage-current lineshapes as a func- 
tion of an applied magnetic field and energy detuning 
(the energy difference between (1,1) and (0,2) charge 
states). In the limit of a strong local spin dephasing 
mechanism or nonuniform magnetic field, we then obtain 
a single simple closed-form analytical expression giving 
a full two-dimensional map of the current as a function 
of detuning and magnetic field. In Sec. IIVI we conclude 
with a summary of the main results and a discussion of 
extensions and possible future work. 
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FIG. 1. (Color online) (a) A typical transport scenario in 
the sequential-tunneling regime. In the high-bias limit (<5/i > 
\E a — E G | , T, where Ej is the energy of level j) electrons hop 
only in one direction: from the left lead to the dot at a rate 
oc Tl, and from the dot to the right lead at a rate oc Fr. 
(b) The transport cycle in the sequential tunneling regime is 
a loop that carries the dot from an iV-electron state (labeled 
with {a, a'}) to an N + 1-electron state (labeled here with 
Greek letters {a,/3}). A current blockade will be set up if 
one or some of the accessible N + 1-electron states a have 
negligible transition rates W aa — 0, in which case the current 
may be limited by rates Wp a inducing transitions within the 
N + 1-electron subspace. 



II. LEAKAGE CURRENT THROUGH 
BLOCKADED SYSTEMS 

In this section we establish a generic theory of the cur- 
rent through blockaded nanostructures, incorporating in- 
ternal transition rates between levels that can lead to a 
lifting of the blockade. The formalism developed in this 
section will be used in later sections to derive simple ana- 
lytical formulas for the dependence of the current on vari- 
ous parameters in the specific case of the Pauli spin block- 
ade of transport through a double quantum dot, account- 
ing for inelastic cotunneling processes. Throughout this 
paper, we work in units where <?/zb = e = ks = h = 1, 
with electron g- factor g, Bohr magneton /xb, electron 
charge e, Boltzmann's constant fee and Planck's constant 
H. 

We consider a nanoscale system weakly tunnel coupled 
to leads, set at chemical potentials /i; and with tunnel- 
ing rates T/ (I = L, R for the tunneling rate between 
the system and the left and right leads, respectively, see 
Fig. HJa)). In the sequential-tunneling regime, electron 
charge is transported from left to right through energy- 
conserving transitions between TV-electron states of the 
system, denoted a, and N + 1-electron states, denoted 
a. These transitions are associated with the exchange 
of an electron with one of the leads (Fig. [Tfb)). Fi- 
nally, we consider the simplifying limit of a large bias 
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Sfx = hl — p R > 0: 

I Ml - /w| > T, (1) 

Imjj - /Vol > r, (2) 

with chemical potential /i aff = — E a (where is the 
energy of isolated system level i and Latin characters are 
taken to run over all system eigenstates, independent of 
the occupation number, i.e.: i = {a, a}). In this limit, 
energy-conserving transitions that add one electron to 
the system a — > a (with rate W a<7 ) necessarily involve the 
removal of an electron from the left lead and transitions 
that remove an electron from the system, a — > a (with 
associated rate W aa ) , involve the addition of an electron 
to the right lead (see Fig. QJa)). Other processes that 
change the electron number are exponentially suppressed. 

The sequential-tunneling current in the high-bias limit 
(defined by Eqs. (fTJ) and ^) is given simply by 

aa 

where pi solves the (stationary) Pauli master equation 
for the diagonal elements of the system density matrix: 

^=J2 W ijPo-WiPi = 0; W i = Y t W iii (4) 

3 3 

where we define the average p(t) = lim T _ ) . 00 i f Q T dtp(t). 
The diagonal elements of the stationary density matrix 
pi must satisfy the normalization 

i 

Use of the classical (Pauli) master equation to describe 
the diagonal elements of the system density matrix is 
strictly valid in the high-bias, weak-coupling limit \Sp\ > 
Fl,Ri where coherences (off-diagonal elements with re- 
spect to the isolated system energy eigenbasis) decay to 
zero on a time scale ~ 1/1 I th a * i s short compared to 
the tunneling time ~ 1/T;. In Eq. ([3]), we have explicitly 
assumed that higher-order current-carrying cotunneling 
corrections cx T^Tr are small relative to the sequential- 
tunneling terms a b. In Sec. IIII1 we will account for 
cotunneling processes that do not carry current, cx Tf, 
involving exchange with the same lead. These processes 
typically dominate over the current-carrying cotunneling 
processes in the case of a double quantum dot considered 
in Sec. UTT] 

Solving the linear system given by Eqs. (jH) and ([5]) for 
the stationary populations p a is sufficient to determine 
the current / from the set of all rates Wij. However, it is 
physically intuitive to switch to new variables ki, defined 
in terms of the current J, stationary populations pi and 
the total escape rate from state i, Wi = J2j as: 

Wipi flux out of state i 
I flux into all states a 



The quantities ki have a natural physical interpretation: 
ki is the number of times state i will be visited, on av- 
erage, per transport cycle (a —> a — s- a' —> . . ., depicted 
in Fig. HJb)). From the definition ^ and normalization 
(0, we find the current 

This formula can be understood directly in terms of the 
interpretation given above for the coefficients k,. The 
average time to leave state i if it were occupied is 1/Wi, 
while ki is the number of times (on average) that state 
i is occupied in each transport cycle. The ratio ki/Wi is 
therefore the average time spent in state i per transport 
cycle and so the total average time per transport cycle 
(average time to transfer an elementary electron charge) 
is simply ^2^ki/Wi. The inverse of this time is the rate 
at which charge is transfered from left to right, giving 
the current, Eq. (J7J). 

The coefficients ki can be found systematically in terms 
of the rates Wij using the identity 

ki = J2 p ijkj; ii ,,, u - . (8) 
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where the branching ratios, Pij, give the probability for 
a transition to state i conditioned on starting in state 
j. Eq. ||HJ) follows directly from Eq. and the defi- 
nition ki cx Wipi (Eq. ©). Solving the linear system 
given in Eq. ([5]) and substituting the result into Eq. Q 
is formally equivalent to solving for the populations pi 
and substituting the result into Eq. ^ for the current. 
However, for the particular case of blockaded systems, 
we will find that Eq. © lends itself better to approxi- 
mation schemes and often the solution for the ki can be 
determined quickly on physical grounds without directly 
solving the linear system. 

III. INELASTIC COTUNNELING AND THE 
PAULI SPIN BLOCAKDE REGIME 

Here we apply the formalism of Sec. |TT] to perform 
an explicit microscopic calculation for the leakage cur- 
rent through a double quantum dot in the Pauli spin 
blockade regime, accounting for transition rates due to 
inelastic cotunneling processes. Inelastic cotunneling is a 
second-order tunneling process associated with a change 
in energy of the isolated quantum-dot state with a com- 
pensating change in energy of the lead state. Since the 
total energy is conserved, the energy of the combined dot- 
plus-leads system is, of course, unchanged in this process. 

The simple formalism derived in the previous section 
allows us to efficiently obtain closed-form analytical ex- 
pressions for the leakage current in terms of all transi- 
tion rates Wij . We begin by reviewing the calculation of 
transition rates due to sequential tunneling and inelastic 
cotunneling, presented in ref. l25l . 
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FIG. 2. (Color online) Isolated double-quantum-dot en- 
ergy eigenstates and transition rates between them. The 
spin-polarized triplet states \T±) are split from the spin- 
unpolarized triplet |Tb) by an applied magnetic field B and 
the hybridized spin-singlet states \S±) are split from |Tb) by 
the detuning-dependent quantities u±, given in Eq. (|35p . be- 
low. 



Transition rates 



Our starting point is a Hamiltonian, 



H = Hn 



(9) 



where Hq gives the free Hamiltonian of the double quan- 
tum dot and Fermi-liquid leads. In the subspacc of (1, 1) 
and (0, 2) charge states, the isolated double-dot Hamilto- 
nian eigenstates consist of two spin-singlets \S±), describ- 
ing (1, 1) and (0, 2) charge states hybridized due to inter- 
dot tunneling t, as well as three spin triplets (|T±) , |T )), 
see Fig. [2j Sequential transport further involves the 
(0, 1) charge state with spin a =t,|, denoted: \a). The 
Hamiltonian H^i describes the coupling of lead I to dot 
I (I = L, R for the left and right dot and lead, respec- 
tively), with tunneling amplitude U (a full specification 
is given in the appendix, below). 

Transition rates are calculated via Fermi's golden rule: 

W km = 27rY,P(*)\(fk\V\im)\ 2 S(£ fk ~£ m ), (10) 

if 

where i, f label the initial and final states of the leads, 
respectively, m, k label the initial and final states of the 
double dot, p(i) describes a product of initial Fermi dis- 
tributions in the leads, held at electron temperature T 
and chemical potentials /xj, and £i m ,£fk give the total 
initial and final energies of both the double-dot and leads. 
The effective perturbation V , accounting for transitions 
up to second order in H d i, is 



v = J2 H ^+ E i a >< 6 i 

l abgll' 



(qj H dl |g) (gj H AV \b) 
f - f 



[11) 



where the indices a, b, g describe the collective state of 
the double dot and leads. 

In the high-bias regime (Eqs. (QJ and ((2])), we find the 
sequential tunneling rates from Eqs. ITUlandlTTI 



W aa = T L J2 \Al 



o\2 



aa |2 



(12) 
(13) 



where the tunneling rates T/ and transition matrix ele- 
ments A\ 3 are, respectively, 



= 2ttui\U\ 2 , 

4f = 0'l4li'> 



(14) 
(15) 



Here, vi is the density of states per spin at the Fermi 
level in lead I . We assume that both vi and tj are approx- 
imately energy-independent in our regime of interest 4i 
The operator d\ a creates an electron in single-particle or- 
bital I with spin a. 

The inelastic cotunneling rates W a p, arising from the 
second-order term in Eq. (fTT| with I = I', are given by2£ 



W a/3 = 2cTM aP F{ujp a /T), 



(16) 



with dimcnsionless prcfactor c, matrix elements M a p 1 
and energy-dependent factor F(uj/T) given by 



E 



A - 



A + e — 2U' — 26 fx 



(17) 
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/3cr 
La' 



EKi 2 i<'i 2 (i8) 



F(u/T) 



lo/T 



1 



,-u/T ' 



(19) 



The coefficient c ~ Tf in Eq. (|T7|) reflects the second- 
order nature of the cotunneling process and we have 
taken the convention /i£ = 0, [j,r = —5 fx with positive 
bias £/i > 0. Here, the energy detuning e = E(iu — -E(o,2) 
measures the separation in energy between (1,1) and 
(0,2) charge configurations, and A = .En i) — U sets 
the energy of the (1,1) charge configuration with on-site 
charging energy U. The transition matrix elements M a p 
arise from processes involving spin exchange with the left 
lead (associated with virtual states in the (0, l)-subspace) 
or spin exchange with the right lead (associated with 
virtual states \j) in the (1, 2)-subspace). The dominant 
energy and temperature dependence of the cotunneling 
rates is due to the function F(lj/T), which arises from an 
integral over Fermi functions fi(E) = l/{eS E ~^/ T + 1): 



F(w/T)=i ^ dEfi(E)[l-ME- 



>)} 



(20) 



In our analysis, we have neglected resonant cotunneling 
contributions, which formally lead to a divergence in eval- 
uating rates directly from Eq. (fTT| . However, these con- 
tributions can be systematically regularized^ and are 
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suppressed exponentially in the high-bias limit consid- 
ered here. We have further neglected current-carrying co- 
tunneling processes (those arising from the second-order 
term with I ^ I' in Eq. (fTT]) h We find that these pro- 
cesses are suppressed relative to the considered processes 
by at least a factor ~ U'/U, where U' is the nearest- 
neighbor charging energy^ 

Since the energy dependence of the rates W a p will play 
an important role in the following analysis, it is useful 
to consider F(u>/T) in the limits of large positive and 
negative energy difference ui at low T: 



F(u/T) ~ -Biu/T), 



F(u/T) ~ |w/T|e 



-\ui/T\ 



uj > T, 
uj < -T. 



(21) 
(22) 



Eq. (|2ip reflects the fact that the inelastic relaxation 
rates increase for large energy-level separation uj > T, as 
the density of states of the environment increases and Eq. 
(|22[) describes exponential suppression of excitation com- 
pared to relaxation processes, consistent with detailed 
balance. 

In the high-temperature limit, the inelastic cotunncl- 
ing rates for both excitation and relaxation approach a 
constant, energy-independent value since, in this limit, 



F(w/T) ~ 1, T > M 



(23) 



This high-T limit has been explored in the context of 
Pauli spin blockade in previous works i 17 ' 24 In the present 
work, we are more concerned with the limits where the 
energy dependence described by Eqs. ([21]) and (|22|) is 
significant in determining the leakage current^ 



B. Leakage current: No local dephasing 

Solving the linear system (Eq. ([H]l) for fc, with the rates 
given in Eqs. (fT2j). (fl"3"|) . and ([To]) immediately gives the 
current via the expression in Eq. (J7]) . For any set of pa- 
rameters, one can find the leakage current by solving the 
full linear system, giving a complex expression in general. 
However, to understand the physical significance of the 
results, or in order to perform experimental fits to traces 
of leakage current vs. magnetic field or energy detuning, 
it is useful to derive simple analytical expressions, valid 
in experimentally relevant limits. In this section, we de- 
rive expressions for the leakage current in the limit where 
there is no significant local dephasing mechanism, lead- 
ing to decay rates that are comparable for all three spin- 
triplet states, Wt + ~ Wt_ ~ Wt - This limit applies to 
double dots in silicon or carbon when there is no mag- 
netic field gradient. In Sec. IIII CI below, we consider the 
opposite limit of a strong local spin dephasing mechanism 
or nonuniform magnetic field, leading to Wt 3> Wt± , in 
which case the analysis simplifies considerably. 




FIG. 3. (Color online) Current I vs. magnetic field B at zero 
detuning (e = 0) from the full expression given in Eq. I|28[l 
(black solid line) and the approximation, Eq. (|30|) (red cir- 
cles). This form is valid when the temperature T is large com- 
pared to the tunnel coupling t and when the direct dot-lead 
tunneling rates T; are large compared to inelastic cotunneling 
rates W a p, i.e. T > t and F > W a p. 



1. B-field dependence (high-T limit) 

First we restrict ourselves to the dependence of the 
leakage current on an applied magnetic field B, which 
splits the spin-polarized triplet states \T±) from \Tq) (see 
Fig. [2|). We further consider the limit T; 3> W a p at zero 
detuning, e = 0, in which the current is dominated by 
rates of escape from the three spin-triplet states Wr a *C 
Ws± , W a . Noting that ks ± ,k a ~ 0(1) in this limit since 
the singlets \S±) and one-electron states \a) are accessed 
at most once per transport cycle, Eq. (0 simplifies to 



( k 



To 



+ 



kT_ 



\W To W T+ W T _ 



T t » W a p. (24) 



We obtain the relevant coefficients k a from Eq. @, 
with the cotunneling rates given in Eq. ([T!)|) . We find 
k a = 1/2, since both spin states, a =t,i, arc equally 
probable and one of the two is accessed in each trans- 
port cycle. For an unpolarizcd source lead, we find the 
branching ratios -Pr+t = -Pr_| = 1/2, Pt c = 1/4, and 
finally, directly inserting the relevant rates in the limit 
e = 0, t < T, we find Pt t± =1/2. Inserting these re- 
sults into Eq. ((5J) and solving the linear system in terms 
of the two remaining branching ratios gives: 

k T 



1 

4 



2 - P+ - P- 
P± 

1-P+-P- 



with branching ratios 



P± 



T±T 



W T± T 



Th- 



ro 



(25) 
(26) 

(27) 



Substituting the cotunneling rates given in Eq. (|16[) into 
Eqs. (|27p. (|23|) . and (|2^|) . and inserting the results for 
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e/T 



FIG. 4. (Color online) Detuning dependence at B — from 
Eq. (011) for a = 100 and t = 5T (black solid line), t = 3T 
(blue dashed line), and t = T (red dotted line). The central 
peak with width ~ t is described by the function g(e,t,T) 
and the broad background describes a Lorentzian of width 
Se — at. 



T > t into Eq. (|24jl directly gives an expression for the 
current vs. B: 



with 



G{x) = 3 



2(cosha; — 1) + xsinhx 
2(cosh:r — 1) + 5xsinhx 



(28) 



(29) 



Since G(B/T) differs from a constant only at third order 
in B/T, while I(B 7 e = 0) is exponentially suppressed 
for B/T > 1, to a very good approximation we take 
G(B/T) ~ G(0) = 1 leaving the simple expression 



x 4 BIT 



T>t,r l ^>w a 



. (30) 

The approximate expression, Eq. (|30[) , is virtually indis- 
tinguishable from the full expression given in Eq. (|28[). 
see Fig. O 

At zero magnetic held, B = 0, Eq. ([50)1 simply gives 
/ ~ nWcot, where n = 4/3 gives the average number 
of electrons that pass through the double dot between 
"blocking events". This number is 4/3 if three of four 
(1,1) charge states block current - in this case, the three 
spin triplets - see ref. H^. W co t = C T is the rate at which 
any one of the triplets converts to a singlet through a 
cotunncling process in the limit T > t. 

Physically, the current in Fig. [3] falls to zero when 
\B\ > T since the excitation rate out of the ground-state 
triplet is exponentially suppressed. 



2. Detuning dependence 

To find the full e-dependence of the current, even in the 
limit Ti ^> W a p we find that it is necessary to include 



the escape rates from the singlets Ws± , which can control 
the resonant-tunneling current at large detuning, where 

w s± ~w Toi w T± . 

As in the last section, we aim to find a good approx- 
imate solution for the current, starting from the linear 
equations given in Eq. (JSJ). In the limit T/ 3> W a p, 
assumed here, the escape rate from the singlets can be 
taken to be essentially instantaneous except at suffi- 
ciently large detuning |e| 3> t, since W a s ± « {t/ e ) 2 ^R 
for e — > ±oo (the rate W a s ± is limited by the overlap of 
|5±) with the (0, 2) charge state). It is therefore sufficient 
to approximate rates by their large-e forms to determine 
ks± ■ We assume that transition rates between the sin- 
glets are small compared to the direct escape rate (i.e., 
W S _ S+ - (t/e) 2 c\e\ « mm(W aS+ ,W aS -) « (t/e) 2 T R: 
which simplifies to c\e\ -c r#). Provided this is satis- 
fied, in the limit |e| 3> t, we find that only fcs + (fcs_) 
is relevant for e > (e < 0), allowing us to introduce a 
single parameter ks = 6(e)ks + + 6(—e)ks + , where 9{x) 
is a Heaviside step function. At B = 0, we further find 
that kx + = kx_ = fcr- The remaining three independent 
parameters are then given by the linear equations, from 
Eq. ([5]) after inserting all rates: 



kx 
ks 



1 



1 

4 + 2 + g 
1 2 
4 + l + g 
1 1, 

T + 77 fc Tn 



kr 
kx 



1 



3 + r/ 
1 



3 + ?/ 



ks, 



(31) 
(32) 
(33) 



with 



g = 9{e,t,T) = , (34) 



Ve 2 + Si 2 



Ve 2 + 8t 2 T e 



(35) 



and 



aS 



2T R 
cT 
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T 



lei » t. 



(36) 

Here we have introduced a new energy scale Se, giving the 
value of the detuning at which the inelastic cotunncling 
rates are comparable to the escape rate from the double 
dot: 



Se = at; a = \J 3Tji/8cT. 



(37) 



As the detuning e — > ±oo is increased, the escape 
rate from the relevant singlet state W a s± oc (t/e) 2 Tn 
decreases until it becomes smaller than the cotunncling 
rates ~ cT. In the extreme limit, rj = and g = 1, 
Eqs. (|3"T|) . (|32|) . and (|33|) only have the singular solution 
kx = kx = ks — oo, reflecting the fact that each state is 
visited an infinite number of times; the system becomes 
'stuck' in loops, as depicted between states \a) and \(3) in 
Fig. [U To arrive at the leading finite corrections in the 
limit of large detuning e 3> Se (equivalently, rj <C 1), we 
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FIG. 5. (Color online) Detuning dependence at B — from 
Eq. (gHJ) for t = 5T and a = 200 (black solid line), a = 50 
(blue dashed line), and a = 10 (red dotted line). 



set .9 = 1, and determine the leading asymptotic solution 
of Eqs. dSTJ, (|32J}, and §5$ for small 77, giving 



k T 



fc T+ 



2 
V 

const. 



const., e ^> Se, (38) 
e > Se. (39) 



In the same large- detuning limit, we have the total decay 
rates (at any finite temperature T) 



W T+ 



W T _ 

w s± 



2cT, 
3cT, 



e — > ±oo, 
e — > ±oo. 



(40) 
(41) 



In the opposite limit of small detuning (e <C Se, or 
equivalently r\ ^> 1), Eqs. (|3"Tj) and (|3"2"j) decouple from 
Eq. dHJ). For |e| < i and low temperature T < t, it 
is necessary to keep the ^-dependence. The resulting 
solutions in this limit are 



k T 



1 

4 



1 

2fJ 



k T+ = 



e < 5e 
1 1 



4 ' 4 ff ' 

with corresponding rates given by 

W T± =cT{l+g), 
W To = cT (2 + .g) . 

Combining the above results gives 

^ fcj ( 4 1 3 
cT 



e < <5e, 



^W n ~ dfv ~ Adf \Se) 

o 



e > Se, 



(42) 
(43) 



(44) 
(45) 



(46) 
(47) 



Since the result in Eq. (|46|) vanishes for e <C Se, but 
dominates over Eq. (|4"7| for e ^> £e, we can simply add 
the two results to find the appropriate denominator for 
the current, Eq. ([7]). giving an expression that closely 



0.5 
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FIG. 6. (Color online) Detuning dependence at B — from 
Eq. QHJl for a = 100 and t = 0.5T (black solid line), t = 0.1T 
(blue dashed line), and t = 0.01T (red dotted line). Open 
circles give the equivalent curves from Eq. (|49|l . valid in the 
limit t < T. 



approximates the current everywhere except possibly in 
a small region around e ~ Se. 

The resulting lineshape for the leakage current as a 
function of detuning e is 



I(B = 0,e) 



cT 



g- l (e,t,T) + (e/SeY 



r ; »VF Q(3 . (48) 



In general, the leakage current lineshape may be dom- 
inated by the function g(e,t,T), due to escape from 
the triplets at small detuning e < t, and by a broad 
Lorentzian with width Se, limited by escape from the sin- 
glets at large detuning, e > t. Eq. (|48|) is plotted in Fig. 
[4] for various values of the tunnel coupling t, demonstrat- 
ing the crossover from a narrow central peak dominated 
by triplet relaxation to a broad Lorentzian background, 
when each of the states is visited many times before an 
electron escapes the double dot. The balance between 
broad Lorentzian and peaked resonant tunneling can be 
tuned with the ratio of escape rate ~ to cotunneling 
rate ~ cT, controlled by the parameter a. The evolution 
of the current vs. detuning as a is varied is shown in Fig. 

El 

The lineshape given in Eq. (|48[) simplifies considerably 
in the high-temperature limit T > t, in which case the 
function g(e,t,T) ~ 1, leaving a simple Lorentzian: 



I(B = 0,e) 



cT 



1 + (e/Sey 



-, T>t,Ti^W a p. (49) 



Eq. (|49|) is plotted in Fig. [6] and compared with the 
full expression given in Eq. (|48j) in the relevant high- 
temperature limit. 

Eq. (|49|) is consistent with recent experiments on sili- 
con double quantum dots^S 
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e/T 

FIG. 7. (Color online) Map of the leakage current vs. detun- 
ing, e, and magnetic field, B, in the strong-dephasing limit, 
Wt q > Wt± from Eqs. ((50)) and (f5Tj) . A tunnel coupling 
t = 5T was chosen to generate this plot. Cuts vs. B at finite 
e (dotted line) and vs. e at finite B (dashed line) are shown 
in Figs. [8] and [9] below, respectively. 




0.4 - 

0.2 1 



' -20 -10 10 20 

B/T 

FIG. 8. Cut of magnetic field dependence at finite detuning 
e/T = 50 along the dotted vertical line in Fig. [Jj At finite 
positive detuning e, the magnetic-field dependence shows a 
central peak with width set by ~ T due to inelastic escape 
processes involving the triplet To and excited-state singlet S+ , 
followed by a long slow decay at larger B, with a width B ~ e, 
after which the ground state becomes a spin triplet. 



C. Leakage current in the strong-dephasing limit 

An especially simple and ubiquitous limit occurs when 
there is a strong local dephasing process^ or a mag- 
netic field gradient, allowing rapid escape for the spin- 
unpolarized triplet state |T n ) (i.e., Wt — W a T ^ W a p). 
The advantage of this limit is the absence of closed 'loops' 
that complicated the analysis in Sec. IIII Bl 

At sufficiently small detuning e < 5e, the rate of tran- 
sition from one of the singlet states to a triplet is small 
compared to the singlet escape rate, W a s ± <C W a s± ■ 
Consequently, the spin-polarized triplets \T±) are visited 



1.5 - 



1.0 - 



0.5 - 

J 




-400 -200 200 400 

e/T 

FIG. 9. Cut of the detuning dependence at fixed finite mag- 
netic field B/T = 5, along the dashed horizontal line in Fig. 
[Jj At finite magnetic field, the detuning-dependent leakage 
current due to inelastic cotunneling is asymmetric in e, since 
the ground state has spin-singlet character for e > 0, but 
spin-triplet character for e < 0. 

at most once in each transport cycle, giving kx ± = 1/4, 
since each of the four spin states in the (1,1) charge con- 
figuration has equal probability of being occupied during 
a transport cycle. This gives immediately, from Eq. ([7]), 



W T l + Wt} ' 

r, » w aP , w To » w T± , w aS± » w aS± . (50) 

It is important to emphasize the generality of the simple 
expression given in Eq. ([50]) . In particular, this expres- 
sion is valid for arbitrary spin-flip processes leading to 
transitions from the spin-polarized triplet states to the 
singlets or spin-unpolarized triplet To, e.g. Ws + t± 0, 
Wt t± 0, or direct transitions leading to escape from 
the double dot, W a r ± ^ 0. These processes can be medi- 
ated by coupling to nuclear spins, spin-orbit interaction, 
or any other mechanism. The existence of a fast local de- 
phasing mechanism (e.g., coupling to nuclear spins) even 
when spin relaxation rates may be slow is common. Even 
if there is no dephasing mechanism, a magnetic field gra- 
dient across the double dot is sufficient to reach the limit 
W To > W T± H 

When only inelastic cotunneling processes account for 
the rate Wt± , the direct transitions carrying an electron 
out of the double dot vanish W a T± = and the remaining 
contributions can be found directly from Eq. (fl"6| : 

W T± /cT = {W ToT± + W s+T± + W S _T±) /cT 

= F[±~ + ^ 1 \ /. (5i) 

We recall that the energies w± are definied in Eq. ((35)) 
and the functions F(x) are defined by Eq. (|T5)) . Insert- 
ing the rates given in Eq. (|51[) into Eq. (|50p immedi- 
ately gives a complete map of the leakage current, with 
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amplitude determined by the high-temperature cotunncl- 
ing rate ~ cT and all other features determined by only 
three dimcnsionless parameters: the magnetic field, de- 
tuning, and tunnel coupling, scaled by the temperature: 
B/T,e/T,t/T. We show a map of the leakage current 
as a function of B/T and e/T in Fig. [7] for the case of 
t = 5T. 

Taking the limit B = e = 0, T > t, we find that the 
current saturates at a maximum value: 

I(B = 0, e = 0) = 4cT = nW cot , (52) 
Wcot = 2cT, T » t. (53) 

Here, we find the average number of electrons passing 
through the double dot for each cotunneling event is fi = 
2, in contrast to the situation in Sec. IIII B[ where n = 4/3 
for the same parameter values. This is due to the fact 
that now only two out of four of the (1,1) charge states 
are blocked (|T ± )), giving n = 4/2 = 2. 

A cut of 7 vs. B at finite detuning e > is shown 
in Fig. [5J The central peak in this figure has width 
~ T, due to thermally activated escape from the ground 
state triplet through the To state. The broad background 
current falls to zero at B ~ e, where the ground state be- 
comes a spin-triplet (|T_) [\T + )\ for B > [B < 0]). A 
similar effect is shown as a function of e at finite mag- 
netic field B > T in Fig. |H1 Here, for e < the ground 
state becomes a spin triplet and current is suppressed ex- 
ponentially in B/T, whereas for e > the ground state 
is a spin-singlet and relaxation processes can still lead to 
escape with a slow rate ~ 1/e until e becomes very large. 
A detuning asymmetry such as this one is often ascribed 
to phonon-assisted tunneling, but can result (as it does 
for the inelastic cotunneling mechanism considered here) 
from any other mechanism for which excitation is expo- 
nentially suppressed relative to relaxation. 

1. B- field dependence (low-T limit) 




FIG. 10. T = expression for / vs. B (valid for T < t). 



processes are exponentially suppressed and the system 
becomes locked in the ground-state triplet. The dip at 
B = occurs because relaxation processes from T± to To 
vanish when the levels become degenerate, while at small 
finite B, an additional "escape route" is available for the 
highest-energy triplet through |T ). 

From Eq. ([50j) . it is clear that the current will ex- 
perience a dip at B = whenever the rates Wt± are 
reduced at B = 0. This effect becomes especially pro- 
nounced for contributions to Wt± from, e.g., the spin- 
orbit coupling, which must necessarily vanish at B = 
due to time-reversal invariance^ This effect due to spin- 
orbit coupling has been demonstrated in the context of 
the Pauli spin blockade regime using a phenomenologi- 
cal model that preserves time-reversal, but hybridizes the 
triplet and singlet states*^ For a microscopic theory, an 
additional magnetic-field gradient or local spin dephas- 
ing process is likely necessary to arrive at this conclu- 
sion in general, since the spin triplet state \Tq) does not 
hybridize with the spin singlets at leading order in the 
spin-orbit coupling.— 



In the low-temperature limit, T -C t, we use the ap- 
proximation given in Eq. (|2ip to find the relevant escape 
rates from Eq. ([BTjl . For, e.g., B > and e = 0, these 
rates are 

W T+ =cB + ^(V2t + B^Q(V2t + B), (54) 

(55) 



W T _ = C - (V2t - B^j Q(V2t - B). 



Inserting these rates into Eq. (f50|) directly gives the 
low-temperature expression for the current, previously 
reported in ref. Hfl 

(V2t+\B\) 1 ' U ' 

(56) 

Eq. (f5"6")) is plotted in Fig. [TO] The current falls to zero 
at |B| = \flt when the ground-state triplet falls below 
the ground-state singlet \SJ). At larger B, excitation 



IV. CONCLUSIONS 

We have presented a generic and simple procedure for 
calculating leakage current through blockaded structures. 
Using this generic theory, we found simple analytical ex- 
pressions for current lineshapes as a function of an ap- 
plied magnetic field B, energy detuning e, and inter- 
dot tunnel coupling t. These lineshapes fully account 
for inelastic cotunneling in two limits: weak-dephasing 
and strong-dephasing. The results we have found in the 
weak-dephasing limit are consistent with recent experi- 
ments performed on silicon double quantum dots^ and 
may be applicable to carbon-based double dots, which 
are also expected to have weak spin-orbit interactions 
and only weak dephasing. In the strong-dephasing limit, 
we have found a simple expression that gives the full two- 
dimensional map of leakage current vs. B and e in the 
presence of inelastic cotunneling. Finally, we have given 
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a general understanding of zero-field current dips in the 
limit of strong dephasing. 

To simplify the analysis directly associated with the 
Pauli spin blockade in a double quantum dot, we have 
neglected orbital and valley degeneracy, which may be 
relevant in silicon and carbon-based double dots. Ef- 
fects of these degeneracies can, however, be included in 
a systematic and straightforward way using the general 
methodology outlined in Sec. [TTJ We leave the details of 
such an analysis to future study. 



Appendix: Hamiltonian and eigenstates 

In this Appendix we set the precise definition for the 
Hamiltonian and isolated double-dot eigenstates. The 
starting point is a standard tunneling Hamiltonian for a 
double quantum dot coupled to leads 



H — i?dd + £ Hi + Hdl, 



(A.l) 



It is convenient to define new energy variables 

e = V R -V L -U + U', (A.8) 

A = V R + V L - U - U', (A.9) 
where physically, the energy detuning e gives the relative 
energy difference between (1, 1) and (0, 2) charge configu- 
rations and A describes the absolute 'depth' of the (1, 1) 
charge configuration. Diagonalizing Hdd in the space of 
(1,1), (0,1) and (0,2) charge configurations gives the 
cigcncncrgies, assuming a real positive tunnel coupling, 
t > (and defining E (A) = -U - A): 

E (T = -E (A)-^(e-A)+aB/2, (A.10) 

E T± = -E (A)±B, (A.ll) 
E Ta = -E (A), (A.12) 

E S± = -Eo(A) ~l(eT Ve 2 + 8i 2 ) . (A.13) 
The associated eigenstates are 



where Hdd is the Hamiltonian of the double dot, Hi de- 
scribes Fermi liquid lead I, and Hdi gives the tunnel cou- 
pling between lead I and dot I, with / = L(R) for the left 
(right) dot/lead, respectively: 



Hdd 
H 



Hr; + Hn 



H 7 



e lkacl ka Cika, 

ka 



H, 



dl 



h.c. 



(A.2) 
(A.3) 

(A.4) 



Here, cika annihilates an electron in lead I, orbital state 
k with spin a having energy eika- The operator di a an- 
nihilates an electron in dot orbital I with spin a. The 
Coulomb interaction He inter-dot tunneling Hamilto- 
nian ifr and Zeeman term Hz are 



u . 

—mini - 1) 



Vini 



H T = -tY^ (dl a d RlJ +h.c. 



Hz = 2"]C(nJt ~ ni ^> 



U'n L n R: (A.5) 
(A.6) 
(A.7) 



with number operator defined in the usual way, ni = 
J2cr n i<?' n i<? = dfodio-. In the above expressions, U and 
U' describe the on-site and nearest-neighbor charging en- 
ergies, resepectively, in a constant-interaction model, Vi 
gives the local electrostatic potential for dot orbital I, t is 
the inter-dot tunnel coupling, and B is the applied mag- 
netic field (assumed here to be in-plane so that orbital 
effects are negligible). 
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|0) , (A.14) 

(A.15) 
(A.16) 

(A.17) 

\S±) = VC± \S(1, 1)) T VCt \S(0, 2)) . (A.18) 



\T+) = 4 t <4 t |0) , 
\T-)=d{^\0), 

|T ) = ^=(4A+4i4t) io)» 



where the hybridization of the \S(1, 1)) and |5(0, 2)) sin- 
glet states in |5±) is controlled by the parameters 



C± 



Ve 2 + 8i 2 ± e 



2Ve 2 + 8t 2 



(A.19) 



The singlets are defined more precisely in terms of cre- 
ation and annihilation operators by 



|5(i,i)> = ^(4A 

|S(0,2)) = 4 t 4j0>. 



4A)|0>, (A.20) 
(A.21) 
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